Two-year follow-up of the COVID-19 pandemic in Mexico.

Background: After the initial outbreak in China (December 2019), the World Health Organization declared COVID-19 a pandemic on March 11th, 2020. This paper aims to describe the first 2 years of the pandemic in Mexico. Design and methods: This is a population-based longitudinal study. We analyzed data from the national COVID-19 registry to describe the evolution of the pandemic in terms of the number of confirmed cases, hospitalizations, deaths and reported symptoms in relation to health policies and circulating variants. We also carried out logistic regression to investigate the major risk factors for disease severity. Results: From March 2020 to March 2022, the coronavirus disease 2019 (COVID-19) pandemic in Mexico underwent four epidemic waves. Out of 5,702,143 confirmed cases, 680,063 were hospitalized (11.9%), and 324,436 (5.7%) died. Even if there was no difference in susceptibility by gender, males had a higher risk of death (CFP: 7.3 vs. 4.2%) and hospital admission risk (HP: 14.4 vs. 9.5%). Severity increased with age. With respect to younger ages (0-17 years), the 60+ years or older group reached adjusted odds ratios of 9.63 in the case of admission and 53.05 (95% CI: 27.94-118.62) in the case of death. The presence of any comorbidity more than doubled the odds ratio, with hypertension-diabetes as the riskiest combination. While the wave peaks increased over time, the odds ratios for developing severe disease (waves 2, 3, and 4 to wave 1) decreased to 0.15 (95% CI: 0.12-0.18) in the fourth wave. Conclusion: The health policy promoted by the Mexican government decreased hospitalizations and deaths, particularly among older adults with the highest risk of admission and death. Comorbidities augment the risk of developing severe illness, which is shown to rise by double in the Mexican population, particularly for those reported with hypertension-diabetes. Factors such as the decrease in the severity of the SARS-CoV2 variants, changes in symptomatology, and advances in the management of patients, vaccination, and treatments influenced the decrease in mortality and hospitalizations.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico.

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

Metagenomic analysis reveals differences in the co-occurrence and abundance of viral species in SARS-CoV-2 patients with different severity of disease.

BACKGROUND: SARS-CoV-2 infections have a wide spectrum of clinical manifestations whose causes are not completely understood. Some human conditions predispose to severe outcome, like old age or the presence of comorbidities, but many other facets, including coinfections with other viruses, remain poorly characterized. METHODS: In this study, the eukaryotic fraction of the respiratory virome of 120 COVID-19 patients was characterized through whole metagenomic sequencing. RESULTS: Genetic material from respiratory viruses was detected in 25% of all samples, whereas human viruses other than SARS-CoV-2 were found in 80% of them. Samples from hospitalized and deceased patients presented a higher prevalence of different viruses when compared to ambulatory individuals. Small circular DNA viruses from the Anneloviridae (Torque teno midi virus 8, TTV-like mini virus 19 and 26) and Cycloviridae families (Human associated cyclovirus 10), Human betaherpesvirus 6, were found to be significantly more abundant in samples from deceased and hospitalized patients compared to samples from ambulatory individuals. Similarly, Rotavirus A, Measles morbillivirus and Alphapapilomavirus 10 were significantly more prevalent in deceased patients compared to hospitalized and ambulatory individuals. CONCLUSIONS: Results show the suitability of using metagenomics to characterize a broader peripheric virological landscape of the eukaryotic virome in SARS-CoV-2 infected patients with distinct disease outcomes. Identified prevalent viruses in hospitalized and deceased patients may prove important for the targeted exploration of coinfections that may impact prognosis.

Clinical and Virological Features of Patients Hospitalized with Different Types of COVID-19 Vaccination in Mexico City.

Coronavirus disease 2019 (COVID-19) vaccines effectively protect against severe disease and death. However, the impact of the vaccine used, viral variants, and host factors on disease severity remain poorly understood. This work aimed to compare COVID-19 clinical presentations and outcomes in vaccinated and unvaccinated patients in Mexico City. From March to September 2021, clinical, demographic characteristics, and viral variants were obtained from 1014 individuals with a documented SARS-CoV-2 infection. We compared unvaccinated, partially vaccinated, and fully vaccinated patients, stratifying by age groups. We also fitted multivariate statistical models to evaluate the impact of vaccination status, SARS-CoV-2 lineages, vaccine types, and clinical parameters. Most hospitalized patients were unvaccinated. In patients over 61 years old, mortality was significantly higher in unvaccinated compared to fully vaccinated individuals. In patients aged 31 to 60 years, vaccinated patients were more likely to be outpatients (46%) than unvaccinated individuals (6.1%). We found immune disease and age above 61 years old to be risk factors, while full vaccination was found to be the most protective factor against in-hospital death. This study suggests that vaccination is essential to reduce mortality in a comorbid population such as that of Mexico.

Comparative Dynamics of Delta and Omicron SARS-CoV-2 Variants across and between California and Mexico.

Evolutionary analysis using viral sequence data can elucidate the epidemiology of transmission. Using publicly available SARS-CoV-2 sequence and epidemiological data, we developed discrete phylogeographic models to interrogate the emergence and dispersal of the Delta and Omicron variants in 2021 between and across California and Mexico. External introductions of Delta and Omicron in the region peaked in early July (2021-07-10 [95% CI: 2021-04-20, 2021-11-01]) and mid-December (2021-12-15 [95% CI: 2021-11-14, 2022-01-09]), respectively, 3 months and 2 weeks after first detection. These repeated introductions coincided with domestic migration events with no evidence of a unique transmission hub. The spread of Omicron was most consistent with gravity centric patterns within Mexico. While cross-border events accounted for only 5.1% [95% CI: 4.3-6] of all Delta migration events, they accounted for 20.6% [95% CI: 12.4-29] of Omicron movements, paralleling the increase in international travel observed in late 2021. Our investigations of the Delta and Omicron epidemics in the California/Mexico region illustrate the complex interplay and the multiplicity of viral and structural factors that need to be considered to limit viral spread, even as vaccination is reducing disease burden. Understanding viral transmission patterns may help intra-governmental responses to viral epidemics.

Antibody and Memory B-Cell Immunity in a Heterogeneously SARS-CoV-2-Infected and -Vaccinated Population.

Global population immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is accumulating through heterogeneous combinations of infection and vaccination. Vaccine distribution in low- and middle-income countries has been variable and reliant on diverse vaccine platforms. We studied B-cell immunity in Mexico, a middle-income country where five different vaccines have been deployed to populations with high SARS-CoV-2 incidences. Levels of antibodies that bound a stabilized prefusion spike trimer, neutralizing antibody titers, and memory B-cell expansion correlated with each other across vaccine platforms. Nevertheless, the vaccines elicited variable levels of B-cell immunity, and the majority of recipients had undetectable neutralizing activity against the recently emergent omicron variant. SARS-CoV-2 infection, experienced before or after vaccination, potentiated B-cell immune responses and enabled the generation of neutralizing activity against omicron and SARS-CoV for all vaccines in nearly all individuals. These findings suggest that broad population immunity to SARS-CoV-2 will eventually be achieved but by heterogeneous paths. IMPORTANCE The majority of studies on SARS-CoV-2 vaccine-elicited immunity and immune evasion have focused on single vaccines corresponding to those distributed in high-income countries. However, in low- and middle-income countries, vaccine deployment has been far less uniform. It is therefore important to determine the levels of immunity elicited by vaccines that have been deployed globally. Such data should help inform policy. Thus, this paper is very much a "real-world" study that focuses on a middle-income country, Mexico, in which five different vaccines based on mRNA, adenovirus, and inactivated-virus platforms have been extensively deployed, while (as documented in our study) SARS-CoV-2 variants with increasing degrees of immune evasiveness have propagated in the Mexican population, culminating in the recent emergence of B.1.1.529 (omicron).

Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico.

In this study, we analyzed the sequences of SARS-CoV-2 isolates of the Delta variant in Mexico, which has completely replaced other previously circulating variants in the country due to its transmission advantage. Among all the Delta sublineages that were detected, 81.5 % were classified as AY.20, AY.26, and AY.100. According to publicly available data, these only reached a world prevalence of less than 1%, suggesting a possible Mexican origin. The signature mutations of these sublineages are described herein, and phylogenetic analyses and haplotype networks are used to track their spread across the country. Other frequently detected sublineages include AY.3, AY.62, AY.103, and AY.113. Over time, the main sublineages showed different geographical distributions, with AY.20 predominant in Central Mexico, AY.26 in the North, and AY.100 in the Northwest and South/Southeast. This work describes the circulation, from May to November 2021, of the primary sublineages of the Delta variant associated with the third wave of the COVID-19 pandemic in Mexico and highlights the importance of SARS-CoV-2 genomic surveillance for the timely identification of emerging variants that may impact public health.

Genomic Characterization of SARS-CoV-2 Isolated from Patients with Distinct Disease Outcomes in Mexico.

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has shown a wide spectrum of clinical manifestations ranging from asymptomatic infections to severe disease and death. Pre-existing medical conditions and age have been mainly linked to the development of severe disease; however, the potential association of viral genetic characteristics with different clinical conditions remains unclear. SARS-CoV-2 variants with increased transmissibility were detected early in the pandemics, and several variants with potential relevance for public health are currently circulating around the world. In this study, we characterized 57 complete SARS-CoV-2 genomes during the exponential growth phase of the early epidemiological curve in Mexico, in April 2020. Patients were categorized under distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors, the patients were less than 60 years old and with no diagnosed comorbidities A trait-association phylogenomic approach was used to explore genotype-phenotype associations, represented by the co-occurrence of mutations, disease severity outcome categories, and clusters of Mexican sequences. Phylogenetic results revealed a higher genomic diversity compared to the initial viruses detected during the early stage of the local epidemic. We identified a total of 90 single nucleotide variants compared to the Wuhan-Hu-1 genome, including 54 nonsynonymous mutations. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors. IMPORTANCE The genetic association of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with different clinical conditions remains unclear and needs further investigation. In this study, we characterized 57 complete SARS-CoV-2 genomes from patients in Mexico with distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors the patients were less than 60 years old and with no diagnosed comorbidities. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

Dysbiosis and structural disruption of the respiratory microbiota in COVID-19 patients with severe and fatal outcomes.

The COVID-19 outbreak has caused over three million deaths worldwide. Understanding the pathology of the disease and the factors that drive severe and fatal clinical outcomes is of special relevance. Studying the role of the respiratory microbiota in COVID-19 is especially important as the respiratory microbiota is known to interact with the host immune system, contributing to clinical outcomes in chronic and acute respiratory diseases. Here, we characterized the microbiota in the respiratory tract of patients with mild, severe, or fatal COVID-19, and compared it to healthy controls and patients with non-COVID-19-pneumonia. We comparatively studied the microbial composition, diversity, and microbiota structure between the study groups and correlated the results with clinical data. We found differences in the microbial composition for COVID-19 patients, healthy controls, and non-COVID-19 pneumonia controls. In particular, we detected a high number of potentially opportunistic pathogens associated with severe and fatal levels of the disease. Also, we found higher levels of dysbiosis in the respiratory microbiota of patients with COVID-19 compared to the healthy controls. In addition, we detected differences in diversity structure between the microbiota of patients with mild, severe, and fatal COVID-19, as well as the presence of specific bacteria that correlated with clinical variables associated with increased risk of mortality. In summary, our results demonstrate that increased dysbiosis of the respiratory tract microbiota in patients with COVID-19 along with a continuous loss of microbial complexity structure found in mild to fatal COVID-19 cases may potentially alter clinical outcomes in patients. Taken together, our findings identify the respiratory microbiota as a factor potentially associated with the severity of COVID-19.

Genetic Analysis of SARS-CoV-2 Variants in Mexico during the First Year of the COVID-19 Pandemic.

During the first year of the SARS-CoV-2 pandemic in Mexico, more than two million people were infected. In this study, we analyzed full genome sequences from 27 February 2020 to 28 February 2021 to characterize the geographical and temporal distribution of SARS-CoV-2 lineages and identify the most common circulating lineages during this period. We defined six different geographical regions with particular dynamics of lineage circulation. The Northeast and Northwest regions were the ones that exhibited the highest lineage diversity, while the Central south and South/Southeast regions presented less diversity with predominance of a certain lineage. Additionally, by late February 2021, lineage B.1.1.519 represented more than 89% of all circulating lineages in the country.

Emergence and spread of the potential variant of interest (VOI) B.1.1.519 of SARS-CoV-2 predominantly present in Mexico.

SARS-CoV-2 variants emerged in late 2020, and at least three variants of concern (B.1.1.7, B.1.351, and P1) have been reported by WHO. These variants have several substitutions in the spike protein that affect receptor binding; they exhibit increased transmissibility and may be associated with reduced vaccine effectiveness. In the present work, we report the identification of a potential variant of interest, harboring the mutations T478K, P681H, and T732A in the spike protein, within the newly named lineage B.1.1.519, that rapidly outcompeted the preexisting variants in Mexico and has been the dominant virus in the country during the first trimester of 2021.

Genomic Analysis of Early SARS-CoV-2 Variants Introduced in Mexico.

The coronavirus disease 2019 (COVID-19) pandemic has affected most countries in the world. Studying the evolution and transmission patterns in different countries is crucial to enabling implementation of effective strategies for disease control and prevention. In this work, we present the full genome sequence for 17 SARS-CoV-2 isolates corresponding to the earliest sampled cases in Mexico. Global and local phylogenomics, coupled with mutational analysis, consistently revealed that these viral sequences are distributed within 2 known lineages, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage A/G, containing mostly sequences from North America, and lineage B/S, containing mainly sequences from Europe. Based on the exposure history of the cases and on the phylogenomic analysis, we characterized 14 independent introduction events. Additionally, three cases with no travel history were identified. We found evidence that two of these cases represented local transmission cases occurring in Mexico during mid-March 2020, denoting the earliest events described for the country. Within this local transmission cluster, we also identified an H49Y amino acid change in the Spike protein. This mutation represents a homoplasy occurring independently through time and space and may function as a molecular marker to follow any further spread of these viral variants throughout the country. Our results provide a general picture of the SARS-CoV-2 variants introduced at the beginning of the outbreak in Mexico, setting the foundation for future surveillance efforts.IMPORTANCE Understanding the introduction, spread, and establishment of SARS-CoV-2 within distinct human populations as well as the evolution of the pandemics is crucial to implement effective control strategies. In this work, we report that the initial virus strains introduced in Mexico came from Europe and the United States and that the virus was circulating locally in the country as early as mid-March. We also found evidence for early local transmission of strains with a H49Y mutation in the Spike protein, which could be further used as a molecular marker to follow viral spread within the country and the region.